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A LEPTONIC-HADRONIC MODEL FOR THE AFTERGLOW OF GAMMA-RAY BURST 090510 

SOEBUR RAZZAQUE^'^ 
Draft version October 6, 2010 

ABSTRAGT 

We model multiwavelength afterglow data from the short Gamma-Ray Burst (GRB) 090510 using 
a combined leptonic-hadronic model of synchrotron radiation from an adiabatic blast wave. High 
energy, > 100 MeV, emission in our model is dominated by proton-synchrotron radiation, while 
electron-synchrotron radiation dominates in the X ray and ultraviolet wavelengths. The coUimation- 
corrected GRB energy, depending on the jet-break time, in this model could be as low as 3 x 10^^ erg 
but two orders of magnitude larger than the absolute 7 ray energy. We also calculated the opacities 
for electron-positron pair production by 7 rays and found that TeV 7 rays from proton-synchrotron 
radiation can escape the blast wave at early time, and their detection can provide evidence of a 
hadronic emission component dominating at high energies. 

Subject headings: gamma-ray burst: individual (GRB090510) - relativistic processes - shock waves 



1. INTRODUCTION 

] Gamma-Ray Burst science has entered a new era with 
launch of the Fermi 7-ray space telescope. The main 
instrument. Large Area Telescope (LAT), is more sen- 
sitive than any previous instrument in the 20 MeV - 
300 GeV range (Atwood et al. 2009), whereas the Fermi 
Gamma-ray Burst Monitor (GBM) observes the whole 
unocculted sky in the 8 keV - 40 MeV range (Meegan 

', et al. 2009). With the advent of the Burst Alert Tele- 

• scope (BAT), X-Ray Telescope (XRT) and UV-Optical 
Telescope (UVOT) onboard the Swift satellite (Gehrels 
et al. 2004) it is now possible to obtain simultaneous 
multiwavelength data in the optical to multi-GeV 7-ray 
energy range from GRBs. 

GRB 090510 is the first GRB to provide data from si- 
multaneous observtions by Fermi (Abdo et al. 2010) and 
Swift (Hoversten et al. 2009), as well as by a couple of 

; other sateUites. At a redshift z = 0.903 ± 0.003 (Rau 
et al. 2009) the isotropic-equivalent 7-ray energy release 

■ from this short GRB (T90 < 2 s) is E^;,^^ = (1.08±0.06) x 
! 10^^ erg with a fluence of (5.03 ± 0.25) x 10^^ erg cm^^ 

■ in the 10 keV - 30 GeV range (Abdo et al. 2010). Fermi 
LAT detected long-lived emission up to ^ 200 s after trig- 

■ ger (To = 00:23:00 UT, 2009 May 10) in the > 100 MeV 
range. While such high-energy emission, that is tempo- 

' rally extended beyond the keV - MeV emission, was first 
detected in GRB 940217 by the Compton Gamma- Ray 
Observatory (Hurley et al. 1994), this feature is common 
to most GRBs detected with Fermi LAT. Swift XRT 
and UVOT collected data from GRB 090510 between 
~ To -I- 97 s and Tq -I- 1.9 ks before an Earth Occultation 
(EO), and again after Tq -1-5.1 ks (Hoversten et al. 2009). 
Swift BAT collected most data within T90 = 0.3 ± 0.1 s 
(15 - 350 keV), and sparsely between Tq + 0.4 s and 
~ To -h 100 s (Hoversten et al. 2009). 

Smooth temporal evolution of the fiux, F oa i^", of the 
long-lived emission in Fermi LAT (a^ — 1.38±0.07), and 
Swift XRT {ax,i = 0.74 ± 0.03 before EO) and UVOT 
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(«0,i = "0-50^0^13 before EO) observations strongly 
suggest an afterglow origin (De Pasquale et al. 2010). 
Synchrotron radiation by shock-accelerated electrons in 
a decelerating GRB blast wave (Meszaros & Rees 1997; 
Sari et al. 1998) have successfully explained much of the 
broadband afterglow data at radio, optical, and X-ray 
frequencies in the pie-Fermi era. However fitting com- 
bined Fermi and Swift data from GRB 090510 with sim- 
ple e-synchrotron model results in unusual parameter 
values, and most importantly it is difficult to reconcile 
the Fi, oc t~°'v~^ temporal relations (De Pasquale et al. 
2010). More complex scenarios have been proposed to 
model GRB 090510 data such as a radiative fireball in an 
pair dominated environment (Ghirlanda et al. 2009; 
Ghisellini et al. 2009), adiabatic fireball in a low den- 
sity medium and small magnetic field (Kumar & Barniol 
Duran 2009a, b; Gao et al. 2009), and two component 
jet (Corsi et al. 2009). 

Here we present a combined leptonic- and hadronic- 
afterglow model to fit multiwavelength data from GRB 
090510. Inclusion of ion acceleration and radiation in the 
GRB blast wave is a natural and simple extension of the 
e-synchrotron blast wave model, and has been discussed 
by a number of authors (Bottcher & Dermer 1998; Totani 
1998a; Zhang & Meszaros 2001; Wang et al. 2009; Raz- 
zaque et al. 2010). We show that LAT emission in the 
> 100 MeV range is dominated by synchrotron radiation 
from protons accelerated in the external forward shock 
of a decelerating blast wave. (Note that Razzaque et al. 
(2010) considered proton-synchrotron radiation from a 
coasting blast wave.) The XRT and UVOT light curves 
can be reasonably reproduced by synchrotron radiation 
from electrons accelerated in the same external forward 
shock. We present the e- and ion- synchrotron afterglow 
model in Sec. 2, compare this model with GRB 090510 
afterglow data in Sec. 3, and discuss our results in Sec. 
4. 

2. SYNCHROTRON AFTERGLOW MODEL 

With a coasting bulk Lorentz factor of Fq = lO^^Fs 
and an isotropic-equivalent kinetic energy i^k.iso = 
10^^i?55 erg, the deceleration time scale for an adia- 
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batic blast wave in a medium of uniform density n cm ^ 
is (Blandford & McKee 1976; Sari et al. 1998) 

idee « 1.9 (1 + z){E55/nf/^r-^^^ s. (1) 

Later the bulk Lorentz factor evolves as 

r ^ 763 (1 + zf/\E!i5/ny/X^/\ (2) 

where tg (> tdec) is measured in seconds. At the deceler- 
ation time t = tdcc, r « The radius of the blast 
wave, R = AT^ct, is given by 



i? « 1.4 X 10"(1 + z)-^/\E55ts/ny/^ 



cm. 



(3) 



The jet-break time at which T f^i 9^^ (Sari et al. 1999), 
where = 0.1^_i is the jet opening angle, is given by 



t 



jet 



10^ (l + z)(i?55/n)^/3^!./fs. 



(4) 



The fractions of energy injected in a forward 
shock (Blandford & McKec 1976) that channel into 
electrons^ and into ions can be calculated from their 
shock-accelerated spectra. Wc assiimc an electron injec- 
tion spectrum" n'^{y^) oc j'-^ for y„^ ,^ < 7e < 7Lt,e- 
Here = rie(mp/me)T{t) and ^'sa,t,e mini- 

mum and saturation Lorentz factors, respectively, for 
the electrons. In case of ions, we assume an injec- 
tion spectrum n'^iiA) oc ^'X^^ for T{t) < 7a < im,A 
and n'j^iij^) oc i^^^ for i^ j^ < i'a < iLtA- ^^^^ 
7',j ^ = rjA^it) is a break in the spectrum and 7^^^ ^ is 
the saturation ion Lorentz factor. The fraction of shock 
energy carried by the electrons is 

I. _ 1 1 _ fy \k-2 



'k-2 l-ij'hUe)'' 



and that by the ions is 



£A — S,AriA 



ki-1 



-2 ^ 



(6) 



Here and are the number fractions of electrons and 
ions that are accelerated by the shock, respectively, with 
an equal pre-shock number density n = He = tia- 

A fraction es of the shock energy is assumed to gen- 
erate magnetic field, and the magnetic field behind the 
forward shock is given by (Sari et al. 1998) 

B' ^ 297 (1 + ^)3/«e^/'(E55n3)i/8t;3/8 q 

In the regime of our interc ut ( b :s> Cg, the Compton pa- 
rameter F = [— 1 + v^l + 4ee/eB]/2 and the energy 
loss by the electrons is dominated by synchrotron radia- 
tion (Sari & Esin 2001). 

The saturation Lorentz factor for electrons is cal- 
culated by equating the acceleration time to the syn- 
chrotron cooling time in the B' field (eq. [7]) as 



6.8 X 106 tl^^^ 



(1 + 0)-Vl6<i,l/2gV4^^^^^3)l/16 ■ 

^ Refererring to both electrons and positrons. 

* Comoving frame variables are denoted with primes. 



(8) 



Here (j>~^ is the acceleration efficiency for electrons. 
The cooling Lorentz factor, found by equating the syn- 
chrotron cooling time to the dynamic time t'^y^ =tT/{l+ 
z), is given by 



7c, e 



11.5 {l + z} 



-1/8, 



(9) 



For ions, of atomic mass A and charge Z, the satu- 
ration Lorentz factor is calculated by equating the ac- 
celeration time to the shorter of the dynamic time and 
synchrotron cooling time as 



^sat.A ' 



'Ysat.A ' 



2.2 X \{S^{ZlA)t^l''{E^^nflH\l^ 



(1 + ^)1/40^ 
1.2 X 10l"(A/^3/2)^3/16 



; i < id,A 



-; i > td,A- (10) 



(l + ^)3/160y2^V4(^^^^3)l/16' 

The transition takes place at 

td,A « 1.4 X 10^2(1 + ^^^^^-12^-7^-5 g_ (11) 

Note that the cooling Lorentz factor for ions j'^ a ~ 
{A^ / Z^){'mp / rrieY^'^ ^ can be larger than the saturation 
Lorentz factor (eq. [10]). 

To calculate synchrotron spectra at different epoch and 
light curves at different frequencies arising from a for- 
ward shock, it is sufficient to calculate different spectral 
break frequencies and flux normalization along with their 
time evolution (e.g. Sari et al. 1998; Chevalier & Li 2000; 
Panaitescu & Kumar 2000; Granot & Sari 2002). 

The characteristic synchrotron frequencies for the elec- 
trons with the minimum, saturation and cooling Lorentz 
factors respectively arc given by 

hv^^e « 7.7(1 + z)l/2ry2(e5^,5)l/2i-3/2 Ggy^ 
hVsat.e ~ 180(1 + z)-^/Ve"'(^^55/n)^/'i;'/' GcV, 

/ii^ce « 0.5(1 + z)-i/2e-3/2(£;55„2)-i/2^-i/2 (12) 

A transition from the fast-cooling (yc,e < i^m,e) to slow- 
cooling {vc,e > i'm,e) takes place at 

io,e » 1.5 X 10i°(l + z){eBilefnE^^ s. (13) 

In both the fast- and slow- cooling cases the maximum 
e-synchrotron flux is given by (e.g. Sari et al. 1998) 

« 52 (1 zr^U^Buf'^d^iE^, Jy. (14) 

Here ^28/(10^^ cm) is the luminosity distance. Note that 
the synchrotron self-absorption frequency is in the radio 
band (e.g. Panaitescu & Kumar 2000) and we ignore that 
while modeling optical to > GeV data. 

The synchrotron break frequencies for the ions of min- 
imimi and cooling Lorentz factors can be expressed as 
scaling relations to the corresponding break frequencies 
for electrons as 

ym.,A = Z{riA/rigf{mg/mp)^Vm,e, 
v,^A = {A''/Z^){mplm,fu,,,, (15) 

and for the ions of saturation Lorentz factor (eq. [10]) as 

10 4\nE,^f/^ 



hu, 



sat, A 



(1 + z)y^^\tV^ 



TeV; t < td,A 
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Fig. 1. — Modeling of GRB 090510 light curves with proton- 
synchrotron (solid lines) and electron-synchrotron (dashed lines) 
radiation from an adiabatic blast wave decelerating in a uniform 
density medium. The model light curves are computed at 100 MeV 
(black), 15 keV (green), 1 keV (red) and 3 eV (magenta) to be 
compared with the LAT, BAT, XRT and UVOT data points, re- 
spectively. The onset of the model light curves is at < Tq -|- 0.3 s 
for To }t 2400 for a surounding medium density of n = 3 cm~^. 
The other model parameters are i?k,iso = 2 X 10^^ erg, lO"**, 
ep 0.5, eg = 0.3. See main text for more details. 



iysat,A ^ {A/Z)'^{mp/me)i(l)e/4'A)Vsat,e; t > td,A- (16) 

Note that the ion-synchrotron spectrum is always in the 
slow-cooling regime {i^c.A > Vm,A) as opposed to the e- 
synchrotron spectrum which can be in the fast-cooling 
regime early and changes to the slow-cooling regime later. 
The maximum ion synchrotron flux is 



77 max 



ki - 1 



fci-i 



Va 



1 



(17) 



for fci ^ 1 and k2 > 2. 

3. MODELING GRB 090510 AFTERGLOW DATA 

Figure 1 shows light curves at different energies from 
the combined leptonic-hadronic model of a decelerating 
adiabatic blast wave in constant density medium. With 

< 



an initial Fq > 2400, the blast wave decelerates at 



[1]) for the parameters used here: Ei 



55 



0.3 s (eq. 

n « 3 cm^-^. 

We model LAT emission in the « 100 MeV - 4 GeV 
range in Fig. 1 from p-synchrotron radiation {A = 1) 



for which m 



m,p 



< Vp < V, 



c,p- 



In order to reproduce the 



flux decay index — 1.38 ± 0.07 in this range for 
slow-cooling spectrum, one requires fc2 — (4/3)Q!.y + 1 = 
2.84 ± 0.09. The corresponding spectral index is /3j — 
(fc2-l)/2 = (2/3)a^ = 0.92±0.05. For < Vm,p < Vc,p, 
the p-synchrotron flux scales as F^, oc t^/^v^/"^. A con- 
straint on hvm,p ^ 100 MeV as early as w Tq -|- 0.4 s re- 
quires that •qp<bx 10^ (eq. [15]) with es ~ 0.3. These 
requirements together with the flux level (eq. [17]) needed 
to reproduce LAT data constrain the p-synchrotron ra- 
diation component. The fraction of jet energy needed in 
shock-accelerated protons (eq. [6]) is w 0.5 for fci ^ 
and ip w 10-*. The rise of the > 100 MeV LAT flux at 
T < To + 0.4 s could be consistent with the rise before 
the blast wave enters the self-similar regime (Sari 1997). 
The p-synchrotron flux in the optical to X-ray is much 



below the XRT and UVOT data. 

We reproduce the XRT light curve, averaged over 
0.3 - 10 keV range, with decay index ax i = 0.74 ± 
0.03 before the EO at T w Tq 1.43 ks as from e- 
synchrotron radiation. The required electron index is 
k = (4/3)ax,i + 2/3 = 1.65 ± 0.04 for > Vm,e > Vc,e 
in the fast-cooling case, which is valid for a time T < 
To 2 X 10^ s (eq. [13]). Note that the spectral index 
= fc/2 = (2ax,i + l)/3 = 0.83±0.02, is close to that 
of 0j from p-synchrotron radiation. In order to produce 
hvm,,e S 1 keV at the beginning of XRT observation at 
T r^'To + 100 s, we require rje < 20{me/mp) (eq. [12]). 
Together with parameter ~ 5 x 10""* required to pro- 
duce the XRT flux level, we calculate the fraction of jet 
energy in electrons to be €e « 10~^ (eq. [5]) with = 1. 

Electron-synchrotron flux in the UVOT range is in 
the frequency range Vc.e < Ve < i^m,e and scales as 
Fi, oc t^^/*zy~^/^. Although the observed flux fitted with 
aoA = ~0.50to ll (Dc Pasquale et al. 2010) is different, 
we note that the expected value of a = 1/4 is consistent 
with UVOT data in the T-Tq w 600 s - 1.43 ks interval, 
and with the upper limit at T w Tq + 90 s. 

It is clear that the observed X-ray and UVOT flux 
decay indices ax, 2 = 2.18 ±0.1 and ao,2 — l-13io.io: re- 
spectively, after the EO at T > To-|-5.1 ks (De Pasquale 
et al. 2010) are softer than the e-synchrotron emission. 
If the jet break takes place in between T — Tq w 1.4 ks - 
5.1 ks, then the expected decay index for i/g > Vm,e > J^c,e 
is cx t~'^ which is intermediate between ao,2 and ax, 2 
since k = 1.65. Because of the idealized nature of the 
afterglow model and evolution of the blast wave during 
and after the jet break (e.g. Sari et al. 1999), the ob- 
served flux steepening after the EO could still be due to 
a jet break. For 1.4 ks < t-^^t — Tq < 5.1 ks the jet open- 
ing angle is 0.16 < 6'_i < 0.26 (eq. [4]). If the jet-break 
takes place at T > Tq -I- 100 ks, then O-i > 0.8. 

As shown in Fig. 1 the BAT flux is quite noisy and can 
not be reproduced by either e- or p- synchrotron emis- 
sion. A similar conclusion was drawn by De Pasquale et 
al. (2010) based solely on e-synchrotron afterglow model. 
Sporadic emission in the BAT range could be due to cen- 
tral engine activity, working intermittenly at a much re- 
duced emission level than the initial outburst. 

During the early deceleration phase, the soft photon 
density in the GRB blast wave may be large enough to 
induce 77 — ^ e'^e" pair production and photohadronic 
(pj) interactions by protons, and subsequent cascade for- 
mation. The target photon density can be calculated as 
n'^ie') = 2dl{l + z)F^ / {R^cTe') from the synchrotron 
flux, where e' = hv' = hv{l -\- z)/T. 

We calculate the 77 pair production and p7 pion pro- 
duction opacities from their respective cooling time scales 
and the dynamic time scale for the decelerating blast 
wave model of GRB 090510. The opacities for the 7 rays 
with saturation energies, both from the e- and p- syn- 
chrotron emission, and for the protons with saturation 
energies in Fig. 2. The top panel shows the time depen- 
dence (from right to left) of the opacities at the satu- 
ration energy reached at that time. The bottom panel 
shows the opacities vs. the saturation energies reached 
within the same time interval. Thus the whole time in- 
terval of the top panel is squeezed to fit into each of the 
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Fig. 2. — Opacities for the electron-positron pair production at 
the saturation energies of the electron-synchrotron (dashed Unes) 
and proton-synchrotron (dot-dashed lines) emission. Also shown is 
the opacity for photopion (p-f) interaction (sohd hnes). The hori- 
zontal axis in the top panel corresponds to the opacities at diferent 
time, and in the bottom panel it corresponds to the saturation ener- 
gies reached within the same time interval for each hnes. We used 
the same parameters as used in Fig. 1 to calculate the opacities. 

curves in the bottom panel. 
The saturation energies for the e-synchrotron 7 rays 

scale with time as husat,e ~ 115 4>~^U'^^^ GeV (eq. [12]) 
for the same model parameters used in Fig. 1. For the 
p-synchrotron, the saturation 7-ray energy is hi>sat,p ~ 
4.2 (t)p^t-^/'^ TeV (eq. [16]) for T - To < td,p « 3.7 x 
10^^ (pp s (eq. [11]). Thus the 7-ray saturation energies 
decrease with time while the opacities increase due to a 
flux increase of the target photons. However the opaci- 
ties are small to initiate a substantial e"'"e~ pair cascade 
and accompanying radiation. The same is true for photo- 
pion cascade (see, however, Asano et al. 2009), although 
a small fraction of protons above Ep > 300 EeV {(j)p ~ 1) 
can escape as cosmic rays by converting to neutrons. It 
is interesting to note that the saturation proton energy 

decreases with time as Esat,p ~ 741(^~^t7"^^* EeV be- 
cause of a decreasing bulk Lorentz factor F « 923 ts '^^^ 
(eq. [2]), even though the saturation Lorentz factor in- 
creases with time as j'sat.p (eq- [10]). 

4. DISCUSSION AND CONCLUSIONS 

We have fitted the LAT and XRT light curves from 
p- and e- synchrotron emissions, respectively, before the 
EO from an adiabatic blast wave in a constant density 
medium (Fig. 1). The photon index for > 100 MeV 
emission is 1 -f = 1.92 ± 0.05 for the proton index 
k2 = 2.84 ±0.09, and is consistent within Icr of the mea- 
sured value of the power-law component in the 0.9 s - 



1.0 s time interval (Abdo et al. 2010). At earlier time, 
0.6 s - 0.8 s, where the power-law component is signif- 
icant and harder, the agreement is within 3tT. However 
emission below ~ 20 keV in that time interval, which we 
do not model, may contribute to the hardening of the 
power-law component. Our model is compatible with 
photospheric emission which can dominate in the GEM 
and BAT range (e.g. Razzaque & Meszaros 2006; Gao et 
al. 2009; Toma et al. 2010). 

The isotropic-equivalent 7-ray energy release from 
GRB 090510 is 2-4 orders of magnitude larger than the 
typical range of 10^^ - 10^^ erg for short GRBs, derived 
from pre- Fermi era data (Nakar 2007). Our model re- 
quires even larger, 10^^ erg, of isotropic-equivalent en- 
ergy release in the jet (see, however, Freedman & Wax- 
man 2001; Eichler & Jontof-Hutter 2005). Thus the jet 
of GRB 090510 must be strongly beamed for our model 
to be viable. Depending on the jet opening angle de- 
rived in Sec. 3 from our model fit, the absolute jet en- 
ergy is in between « (3 — 7) x 10^"'^ erg if the jet break 
takes place between 1.4 - 5.1 ks. For a jet break at 

> 100 ks, the absolute energy release is > 6 x 10^^ erg. 
Note that the inferred absolute energy of the Fermi-LAT 
GRB sample is generally high, reaching 10^'^ erg in some 
cases (Cenko et al. 2010). We also note that the micro- 
physical parameter ee ~ 10^"* that we derive for GRB 
090510 using the combined p- and e- synchrotron model, 
is much lower than the typically assumed value for short 
GRBs in the range 10~^ - 10~^ which is derived from 
long GRB data in the pre- Fermi era using e-synchrotron 
model only (Nakar 2007). Finally, a confirm detection 
of large-amplitude and short-time variability, synchro- 
nized to the keV - MeV emission, in LAT data will argue 
against their afterglow origin (see, e.g., Dermer & Mit- 
man 1999) as we modeled here. 

A crucial test for the hadronic emission model in the 
afterglows of GRBs at redshift ;< 0.5, thus avoiding ab- 
sorption in the extragalactic background light (Razzaque 
et al. 2009; Finke et al. 2010), may come from ground- 
based TeV telescopes such as MAGIC, VERITAS, HESS 
and HAWC. An extrapolation of the LAT flux in the 

> 200 GeV range, where e-synchrotron radiation flux is 
negligible and internal 77 opacity is small (Fig. 2), should 
produce similar spectra as in LAT and lower flux from 
p-synchrotron radiation. On the contrary, a synchrotron- 
self-Compton model (e.g. Dermer et al. 2000; Zhang & 
Meszaros 2001) is expected to produce a spectral hard- 
ening and flux increase in the TeV range. 

I thank Fermi LAT and GBM, and Swift team 
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J. D. Finke, N. Gehrels, S. Guiriec, V. Pelassa and 
F. Piron, for providing data and useful comments. This 
work is supported by the Ofhce of Naval Research and 
NASA Fermi Cycle II GI program. 
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